Analytical results for the Morse potential in s-wave ultracold scattering: three 

dimensional vs. one dimensional problem 
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Taking advantage of the known analytic expression of the eigenfunctions and eigenenergies of the 
Morse Hamiltonian, explicit expressions are found for the scattering length a and the effective range 
r e which determine the s-wave scattering of ultracold atoms. The effects on a and r e of considering 
the radial coordinate in the interval [0, oo) or in the extended region (—00, 00) are studied in detail. 
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I. INTRODUCTION 

Consider a gas of particles so dilute that binary colli- 
sions are dominant. The collision between two particles 
with thermal momenta p t h = hkth = ^ /XdB is said to 
be in the ultracold regime if the range of the accompany- 
ing interaction is smaller than the de Broglie wave length 
X dB . Under such conditions, as they collide, the particles 
approach each other more closely than their wavelength, 
and details of the interaction become blurred. That is, 
ultracold collisions are expected to be determined by few 
parameters. 

Ultracold collisions may take place in the degenerate 
limit of a gas of particles for which any particle is per- 
manently within a wavelength from other particles, i. e. 
XdB > n 1 ^ 3 w hh n the density of the gas. The search^ 
generatiorPHSl and further study^l f atomic quantum 
degenerate gases have lead naturally to the analysis of 
collisions between ultracold atoms. 

The purpose of this manuscript is to find analytic 
expressions for the s-wave parameters that describe 
such collisions for a potential that exhibits some of the 
main features of an atom-atom interaction: the Morse 
potential 8 



V(r) = £>((!-. 



-ro)\2 
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where D, ft, tq are positive. This potential is repulsive 
for short distances, exhibits a local minimum with depth 
D, width determined by ft, located at r$ and is slightly 
attractive at long distances. Since its proposal it has 
been extensively used to describe anharmonic features of 
the vibrational spectra of diatomic molecules. 

Simple s-wave analytical solutions for the potential can 
be found for bounoPSl and unbound^ states. The key 
is to use an auxiliary mathematical problem where the 
radial coordinate r, that is physically constrained to the 
interval [0, 00), is allowed to vary in (—00,00). In this 
work we analyze the consequences of using this auxiliary 
problem instead of the original one derived directly from 
the Schrodinger equation. 

In general, for a spherically symmetric potential V(r) 
and elastic collisions, the scattering effects at any relative 
momenta p = hk are contained in the partial wave phase 
shifts Se(k). It can be shown that as k — > 0, the s-wave 



phase shift Si = o(k) can be expanded a d 11 * 12 ! 

1 
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a is known as the scattering length and r e as the ef- 
fective range. For other partial waves 5i(k)/k goes to 
zero as k — > 0; s-wave collisions contribute to the scat- 
tering between bosons and distinguishable particles. As 
a consequence, in those cases, collisions in the ultracold 
regime are expected to be isotropic and characterized by 
the scattering length. 

In this article, expressions for a and r e are obtained for 
the Morse Hamiltonian. We begin in Section [TT] by mak- 
ing a brief revision of the bound and unbound eigenfunc- 
tions of the Morse Hamiltonian that vanish as r — > — 00. 
From those unbound functions, the phase shift <5o(fc) is 
explicitly calculated and the scattering parameters a and 
r e are written in an analytical closed form. In Section 
|III[ we study the bound and unbound eigenfunctions of 
the Morse Hamiltonian with the boundary condition that 
nullify u as r — > 0, which is compatible with a radial co- 
ordinate restricted to the interval [0, 00). In an analogous 
way as for the auxiliary problem, the phase shift 8o(k) 
can be calculated and the scattering parameters are im- 
plicitly found. A comparison between the auxiliary and 
the physical system results is then performed. 



II. RADIAL SOLUTIONS FOR THE 
AUXILIARY PROBLEM 



The Schrodinger equation 
2p,dr 2 



+ E- V(r) 



u{r) = 0. 



(3) 



can be related to the stationary dynamics of a one di- 
mensional collision of two particles with reduced mass 
/i and Hamitonian eigenenergy E, or to a three dimen- 
sional s-wave problem for which the radial wavefunction 
has been written in the form R{r) — u(r)/r. Taking V(r) 
as the Morse potential, Eq. [Tl and introducing the vari- 
ables d = y/2j!D/hft, b = ^lE/hft and z = 2de~ fi{ - r ~ r ^ 
a direct calculation shows that the general solution to 
Eq. (|| is 



u b {z) = e- z/2 C 1 z +lb M{\+ib-d,l + 2ib,z) 
+ e- z/2 C 2 z- ib M(± -ib-d,l- 2ib, z), 



(4) 



where C\ and C 2 are constants to be determined and 

"' (P)nz n 



M{ P ,q,z) = J2 



-^ (q)nn\ ' 



(5) 



n=0 



is Rummer's function^ with (p) n the Pochhammer sym- 
bol. It will be useful to know that 

M(p,g,*) = ^eV'-«[l + 0(|*|- 1 )] (6) 

when the real part of z is positive. 



A. Bound States 

The bound states are determined by having E < 0. In 
this case b = iy // 2fi\E\/h(3 = i\b\, and 

u b (z) = e- z ^C lZ -^M^-\b\-d,l-2\b\,z) 

+ e- z l 2 C 2 z + ^M{\ + \b\-d,l + 2\b\ lZ ). (7) 

Since u b should not diverge when z <C 1 (/3r ^> 1), C\ 
must be zero. We now need to apply a second boundary 
condition which will determine the quantization. Solv- 
ing the 3D radial equation would require demanding 
Ub (z(r)) — when r = 0; however, by applying the 
condition where r — > —00, the wave functions and eigen- 
values take a much simpler form which is analytically 
tractable^. We will analyze the consequences of using 
this method in section [TTTJ When r — > — 00, z — > 00 and 
by using Equation mj it is found that 



u b (z) = C 2 e z / 2 z-^ d 



r(l + 2|6|) 

r(| + |6|-d) 



[l + ofl*!- 1 )] 



(8) 

It is worth noting that Ub(z) grows exponentially as z 
grows unless 1/2 + \b\ — d is a negative integer. Since 
Ub{z) should not grow exponentially in that region we 
define — n = 1/2 + |6| — d, where n is a positive inte- 
ger. This condition determines the quantization of the 
energy levels, b n = \b\ = d — n — 1/2, n € 0, 1, 2, . . . or 
E n = -D + hp^/2D/iJL (n + 1/2) - (/j 2 /3 2 /2^) (n + 1/2) 2 . 
Since b n is always positive then n can only take a finite 
number of values for a given d. This means that the 
Morse potential can only hold a finite number of bound 
states. Since the first argument of M turns out to be 
an integer, the solution can be rewritten in terms of La- 
guerre polynomials as 

u(z)=( /3n!2 ^"l ^ e - z / 2 z^L^ b ^(z) 

n[ >~ \T(2\b n \+n + l)J n [) ' 

with which the bound solutions for I = are fully deter- 
mined. 



B. Unbound States 

The unbound states are determined by having E > 0. 
In this case b = y / 2/i\E\/h(3 — \b\, and 

u b {z) = e- z / 2 C\z +l ^M(^+i\b\-d,l + 2i\b\,z) 

+ e- z / 2 C 2 z- lW M{\ - i\b\ -d,l- 2i\b\, zX-10) 

Again, we apply a boundary condition when r — > —00 
which means that z — > 00 where we require that 
Ub(z) — > 0. Using Equation ([6]) it is found that 



u b (z) 



= e z ' 2 z 



C\ 



1 



r(i + 2»|6|) 



ni 



i\b\ -d) 

0(\z\- v )). 



c 2 



r(l-2z|6|) 

ra-i\b\-d) 



(ii) 



As in the bound case Ub(z) grows exponentially with z 
and we can only play with C\ and C 2 to satisfy the con- 
dition Ub(z(r)) — > as r — > —00. For this we look for the 
relationship between C\ and C 2 that nullifies the factor 
with square parenthesis and find that 



C 2 



r(-2i|6|) 

T{\-i\b\-d) 



\ n-2i\b\) y (12) 



where s means the complex conjugate of s. Therefore we 
define A(b) = r(-2i|6|)/r(| - i\b\ - d), so we can satisfy 
the condition with C\ = CbA(b) and C 2 = CbA(b), where 
Cb is a normalization factor that can depend on b. In this 
manner, the solution has the formic 

u b {z) = 2e- z ' 2 C b 

x fH {A(6)z 4|6| M(l + i\b\ -d,\ + 2i\b\, z)}(,13) 

where 5ft means the real part. 

It is important to analyze the asymptotic behavior of 
the solutions since the scattering phase shift depends on 
this. When r — > 00, z —¥ in such way that Equa- 



tion ( 13 ) simplifies to 

u b (z) -> 2C b 'Si[A(b)z^\ 



(14) 



Writing it in terms of r we get that the asymptotic be- 
havior is given by 

u b {z) « 2C b ^{A(b)(2de firo y W \cos(kr) 

+ 2Cb^{A(b)(2de fSro ) llbl }sm(kr), (15) 

where k = \b\/3 is the asymptotic wave number and 3 
means the imaginary part. 

In absence of a potential, the normalized radial wave 
function has the form Uk{r) = y2/7rsin(A:r). The pres- 
ence of the Morse potential also produces an asymptotic 



sinusoidal solution as seen in Equation (15). However, 



the cosine term results in an s-wave phase shift for the 
auxiliary problem 5q X (k) which satisfies 



tan 6, 



(aux) 



:■)? 



(k) 



{A(k/f3)(2de^>Y k/P } 
%{A(k/(3)(2de^oy k /Py 

On the other hand ^{sl/Q^s} = tan(arg(zs)), so 



(16) 



5 { aux) (k) = &rg [i[A(k/ (3) (2de? r °) 



ik/p 



| - arg A(k/P) - | log(2d) - kr (17) 



modulo 7r. Moreover, 



arg^l (|) = ar g r(-j| 



argrfi - «| 



(18) 



Using an expansion for argT(a; - 
phase shift 



iy) we finally get the 



r(<ma:) 



(fc) = -^(7 + ln(2d)+/3r ) + S, 



where 7 is the Euler-Mascheroni constant and 



* = E 



— arctan — — (- arctan — ? 

pn pn j3 (n— d 



(19) 



(20) 



At this point we have full knowledge of the s-wave scat- 
tering phase shift from which, in principle, we can extract 
all the s-wave scattering information for the Morse poten- 
tial, as we will exemplify when we calculate the scattering 
length and effective range. 

We will now proceed by calculating the normalization 
factor. Using the scattering phase shift, we write the 
asymptotic behavior as 

u fe (r) -> 2C b \A(k/(3)\sm(kr + 8 {k)). (21) 

Following Bethe™the normalization of continuum states 
is defined by their asymptotic behavior in which we re- 
quire that 



Uk(r) 



sin(fcr -I- S (k)). 



Therefore, the normalization factor C b is given by 

1 /|6|sinh2^|6|\ 1/2 



C b 



W 



n 



n) 



-1/2 



-v J n 



(22) 

(23) 
(24) 



where v = d — 1/2. 

To finalize this section we analyze the low energy scat- 
tering behavior of the Morse potential. As the parti- 
cles energy tends to zero, the s-wave scattering ampli- 
tude determined by 8q becomes dominant. For low en- 
ergy, Eq. (pi) defines a, the scattering length and r e , the 



effective range. In order to find expressions for them 
we will first calculate the low energy behavior of So(k) 
and afterwards write the expansion ([2]). By identify- 
ing the coefficients of the expansion wc will find the 
parameters we seek. We begin by using the fact that 
arctan x — x — \ +0(x 5 ) to rewrite the series (20) when 
k//3 < 1 as 



a = -m-d) + i]- p 



i^(2) ( i_ d) + | c(3) (_j +0(fc 5 ),(25) 



where ip^ is the polygamma functional and ip = ip(°\ 
Defining two variables 

r i =\(2 1 + \n{2d)+(3r + iP{\-d)) (26) 



P 



and 







the phase shift is rewritten as 

Aaux) 



\,{¥ (2) i\-d) + m)i (27) 



sr x, (k) = -k V +k 3 ^+o(k 5 ) 



(28) 



On the other hand, using the Maclaurin expansion of 
cot Sq we write, 



kcot5o(k) 



kS (k) 



S (k) 



which yields 



k cot So(k) = h k 

V 



2 f V 
3 



O(5 3 ) 



0{k 6 



(29) 



(30) 



Identifying the terms in the previous expression with the 
ones in Eq. (pi) we find that the scattering length is given 
by 

a = ro + i[27 + ln(2d) + V(|-d)], (31) 

while the effective range is 



2 ^ 2 >(|-rf) + 16C(3) 
3° 3/3 3 a 2 



(32) 



The scattering length and effective range as a func- 
tion of the depth of the potential are illustrated in Fig- 
ure [T] for the case r /3 = 4.15 which could correspond to 
two atoms of 6 Li colliding with an electronic state 3 E+ 
when D = 40meV^. One notices that for d « n + 1/2, 
n = 0,1,2,... the scattering length is not well defined 
since its limiting value from the right would be nega- 
tive and diverging while from left it would be positive 







FIG. 1: (Color online) The scattering length and 

effective range obtained from the auxiliary problem as a 

function of d = \f2jID/Kfi for ro/3 = 4.15. 



and diverging. In the nomenclature of scattering the- 
ory that condition is known as the unitarity limit or the 
zero energy resonance^. For those values of d the Morse 
potential is about to support a new bound state. 

As for the effective range, it is always positive with the 
exception rf<l. This condition is not shared by other 
potentials like the square well which admit positive and 
negative values of r e for extended regions of the potential 
depth. We also observe that the r e resonances are located 
to the left of the a resonances where a becomes zero. 



III. RADIAL SOLUTIONS FOR THE PHYSICAL 

PROBLEM: CONSEQUENCES OF INCLUDING 

r < IN THE AUXILIARY PROBLEM 

An auxiliary mathematical problem was used to find 
the analytical results shown in the previous sections. The 
purpose of this section is to understand better the trade- 
offs of replacing the physical problem by the auxiliary 
one. 



A. Bound States 

First of all, the general methodology, described at the 
beginning of last section, when applied to the auxiliary 
problem yields simple analytical solutions. Nevertheless, 
if we allow r to vary only in the [0,oo) interval and de- 
mand Ub(z(r))\ r= Q — that methodology also yields ana- 
lytical solutions that do not reduce to the simple expres- 
sion |9]). These solutions are 

Uj{z) = Ce- Z ' 2 z + ^M{\ + \b j \-d,l + 2\b j \,z). (33) 

Here \bj\, which determine the eigenenergies, are the pos- 
itive roots of the equation in b given by 



M (| + \b\ - d, 1 + 2\b\, 2de /3r °) = 0. 



(34) 



In Figure [2] we compare the energy eigenvalues that 
result in the physical and auxiliary problems for several 



FIG. 2: (Color online) Comparison between the energy 
values that result in the physical and auxiliary problems 
(AE = E( aux "> — E) as a function of the scaled potential 

depth d — ^2fj,D/hf3 and several values of the product 

ro/3, with /3 the inverse of the potential range and ro the 

equilibrium distance of the potential. 



values of the product ro/3. As noticed before, if d is in 
the interval [n+ 1/2, n + 3/2) for n = 0, 1,2, ..., then pre- 
cisely n+1 bound states are supported for the auxiliary 
Morse problem. For the physical problem, this occurs for 
greater values of d. This effect is more evident for small 
values of d and rof3. For instance, for r /3 ~ 1 the Morse 
potential supports no bound states until d > 0.6. For 
ro(3 ~ 4, the first bound state is found for d > 0.5 which 
is the same value (modulo the limited double precision 
of the computer calculation) found in the real problem. 



B. Unbound States 



For small potential depths, d — > 0, the scattering 



length evaluated using Eq. (31) exhibits a logarithmic 



divergence and the free-particle expresions are not ob- 
tained. In order to verify the physical reliability of this 
property, a must be evaluated considering the physical 
boundary condition ui,(z(r — 0)) = 0. Imposing it, a di- 
rect calculation shows that the radial functions now take 
the form 

u b (z) = 2iC b e- z / 2 % {l(6)z i|6| M (| + i\b\ - d, 1 + 2i\b\,z) } 

(35) 
with 

A(b) = z~ m M (| - i\b\ - d, 1 - 2i\b\,z ) , (36) 

and zq = 2de" r ° . In a similar way as we obtained the 
phase shift before we now get 



So(b) 



arg 



J\b\ 



Mb) 



arg 



M(l-i\b\-d,l-2i\b\,z ). (37) 



Notice that, due to the structure of A(b), the factor Zq 
that appears in the phase shift 5q (which gives rise to 




FIG. 3: (Color online) The scattering length as a 

function of d = y/2fiD/hf3 for ro/3 = 4.15 for the 

physical and auxiliary problems. 



the divergence of a in Eq. (311) is now canceled. From 
Eq. (37) a can be calculated by performing numerically 
the limit k ->• of k cot 5 {k), Eq. ^. In Fig. [3] the 
resulting scattering lengths are illustrated and one can 
see that the divergence when d — > is removed for the 
real problem. As for the effective range, it now becomes 
zero as d — > and for larger values of d the differences 
between the values of r e for the physical and auxiliary 
problem result to be less than one percent for all the 
studied cases. 



IV. CONCLUSIONS 

In this work, analytic expressions have been obtained 
that solve the eigenvalue problem of the Morse Hamil- 
tonian under two different boundary conditions. This 



Hamitonian is widely used to model the s-wave anhar- 
monic vibrations of nuclei in diatomic molecules and sup- 
ports a finite number of bound states. It was shown that 
the eigenvalue of the highest excited bound state derived 
from the boundary condition u(z(r — > — oo)) = dif- 
fers significantly from that derived from the condition 
u(z(r = 0)) = for potentials with a range similar to 
the equilibrium position, /3ro ~ 1, at the unitarity limit, 
i. e. with a potential depth D close to the values that 
yield the possibility for the Hamiltonian to support a 
new bound state. Outside this limit, the difference be- 
tween the energy eigenvalues for the auxiliary and the 
physical boundary condition becomes small. This is con- 
gruent with using the former in the standard analysis of 
molecular vibrations. 

We also derived analytical expressions for the phase 
shift in binary collisions both for the auxiliary and the 
physical problem. From them, the most important pa- 
rameters necessary to describe an ultracold collision, that 
is, the scattering length a and effective range r e were eval- 
uated. A divergence of a predicted for very small poten- 
tial depths d — y/2fj,D/H/3 <C 1 for the auxiliary problem 
was removed by imposing the physical boundary condi- 
tion. This analysis illustrates the fact that, even though 
the scattering length is a property that summarizes the 
asymptotic behavior of a wave function at r —¥ oo, it is 
highly influenced by its behavior at the origin. It is im- 
portant to mention that precisely this observation is the 
basis of the theories that use effective potentials to in- 
corporate scattering effects. Perhaps the best well k nown 
example of the latter is the Gross-Pitaevskii equatio n 1 16 ! 17 ! 
that models an ultracold gas of bosons. 
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